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A  FREE  SHEAR  LAYER  PREDICTION  OF  THE  TURBULENT  WAKE 
OF  A  FLAT  PLATE  TOWED  BELOW  A  FREE  SURFACE 


Introduction 

Recently,  Swean  and  Keramidas  (IJ  have  reported  experimental 
results  for  the  turbulent  wake  of  a  flat  plate  towed  below  a  free 
surface.  We  employ,  here,  the  turbulence  modeled  Navler-Stokes 
equations  to  predict  the  experimental  results.  Prandtl's  free 
shear  layer  approximation  is  used  to  model  the  turbulent  kine¬ 
matic  viscosity. 

The  agreement  between  the  predicted  and  measured  velocity 
profiles  is  found  to  be  excellent. 


Governing  Equations 

We  consider  a  flat  plat  (of  infinite  width)  towed 
horizontally  at  a  velocity  U  and  a  depth  d  beneath  the  free 
surface  of  a  fluid.  Figure  I  is  a  schematic  of  the  situation. 

The  z  coordinate  increases  upwards  with  z  -  0  corresponding  to 
the  undisturbed  free  surface.  The  x  coordinate  increases  down¬ 
stream  of  the  plate  with  x  -  0  corresponding  to  the  trailing  edge 
of  the  plate. 

Far  enough  downstream,  the  plate  Induced  axial  flow  velocity 
u  must  satisfy  [2] 
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together  with  the  boundary  conditions 
"0  at  z  ■  0 


(lb) 


Manuscript  approved  August  28,  1984. 


u  ♦  0 


as  z  *  -  m 


(Ic) 


and  Che  momencun  constrainC  [3] 

0 

pU  /  u  da  *  -D  (Id) 

*  v 

—m 

Here,  v  and  p  denote,  respectively,  the  turbulent  kinematic 
viscosity  and  the  density  of  the  fluid;  while  Dv  represents  the 
viscous  drag  per  unit  width  on  the  plate.  To  model  v^,  we  use 
Prandtl's  free  shear  layer  approximation  [4].  For  the  two- 
dimensional  flow  problem  under  consideration,  this  approximation 
specifics  that 

v  -  kU  (2) 

where  ic  is  a  constant. 


Solution  of  the  Governing  Equations 

The  solution  of  equation  (1)  with  v£  given  by  equation  (2) 
is  a  standard  problem  in  heat  conduction  [3,  Chapter  6].  We  find 
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where  uQ(z)  must  satisfy 
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Here,  H(£)  is  the  Heaviside  unit  step  function  defined  by 


H(C) 
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(4) 


2 


and  uQ(z)  is  an  input  velocity  profile  specified  at  the  axial 
location  x  -  xQ  or 

uq(z)  -  u(xq,z)  (5) 

Because  equation  (Id)  represents  a  diffusion  process,  the  far 
downstream  structure  of  the  wake  is  independent  of  the  details  of 
uQ(z).  Hence,  we  assume  uQ(z)  is  given  by  a  Dirac  delta 
function,  6(C),  at  the  plate  depth  d.  Specifically,  we  set 

uQ(z)  -  -U  Qq6  (  z+d)  (6) 

where 

eo  -  Dv/pu2  (7) 

is  the  total  momentum  thickness.  Equation  (6)  clearly  satisfies 
equation  (3b)  as  required. 

Substituting  equation  (6)  into  equation  (3a)  and  redefining 

<  as 

x  -  Y0q  (8) 

we  obtain 


Evaluation  of  y  and  x 

o 

The  velocity  field  uw(x,z)  behind  a  flat  plate  in  an 
infinite  fluid  is  given  by  1/2  of  equation  (9)  with  d  «  0  or 
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Schlichting  [4]  notes  from  experimental  observations  that,  for 
proper  values  of  Y  and  xQ  ,  equation  (10a)  provides  a  valid 
description  of  the  velocity  field  for 


x/0Q  >  100 


(ID 


In  particular,  along  the  plate  centerline  z  «  0 ,  we  find 
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where  we  have  introduced  the  dimensionless  axial  coordinates 


x  -  x/0  ,  x  ■  x  /0  (12b) 

o  o  o  o 

In  Table  1 ,  the  experimental  results  of  Swean  and  Keramidas 
[1]  and  of  Chevray  and  Kovasznay  £  6 ]  for  a  flat  plate  in  an 
infinite  fluid  are  tabulated.  Also  tabulated  is  the  value  of 
Y ( x— xq )  required  In  equation  (12a)  to  produce  the  experimentally 
observed  centerline  velocity.  From  Table  1,  six  values  of  A x  and 
the  corresponding  six  values  of  yAx  can  be  obtained.  These  are 
given  in  Table  2  as  Is  the  resultant  value  of  y •  Averaging  the 
values  of  Y  in  Table  2  gives  a  best  fit  value  for  Y  as 

Y  -  0.0209  (  13) 


Using  this  value  for  y  »  we  tabulate  in  Table  3  the  calculated 

value  of  x-x  from  Table  1.  Also  given  is  the  value  of  x  . 

o  o 

Averaging  the  values  of  x  in  Table  3  gives  a  best  fit  value 


X  -  X  /0 
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12.3 


(  14) 


Comparison  with  Swean's  Experiments 


In  Figure  2,  the  velocity  profiles  predicted  by  equation 

(10)  are  compared  with  the  results  of  Swean  and  Keramidas  for  a 

flat  plate  towed  in  an  infinite  fluid.  The  comparisons  are  made 

for  their  two  most  downstream  measuring  stations  x/0  *  95.8  and 

o 

x/0o  •  172.1  which  correspond,  respectively,  to  x  -  61.7  cm  and 
x  »  110.8  cm.  In  general,  the  agreement  is  seen  to  be  excellent. 

In  Figure  3,  the  velocity  profiles  predicted  by  equation  (9) 
are  compared  with  the  results  of  Swean  and  Keramidas  for  a  flat 
plate  towed  at  a  depth  d  *  5  cm.  Again,  the  comparisons  are  made 
at  this  two  most  downstream  measuring  stations  x  -  61.7  cm  and 
x  «  110.8  cm.  The  agreement  is  seen  to  be  not  nearly  as  satis¬ 
factory  as  for  the  infinite  fluid  case.  Essentially,  the  model 
substantially  underpredicts  the  observed  spread  of  the  wake  and 
the  observed  free  surface  velocity. 

It  is  apparent  from  Figures  (3a)  and  (3b),  however,  that  the 

measured  velocity  profiles  of  Swean  and  Keramidas  do  not  satisfy 

the  momentum  constraint  of  equation  (Id)  if  their  value  of  0  - 

o 

0.644  for  the  towed  flat  plate  in  an  infinite  fluid  is  used. 

That  is 


-/  ird 


-  ©  >  0.644 

z  o 


Unfortunately,  Swean  and  Keramidas  did  not  measure  directly  the 

viscous  drag  on  the  plate.  Thus,  the  appropriate  value  of  0  for 

o 

the  near  free  surface  towing  condition  can  only  be  estimated  by 
using  their  measured  velocity  profiles. 

For  simplicity,  we  approximate  each  of  the  measured  profiles 
by  a  two  straight  line  fit.  In  Figure  (3a),  we  find  the  three 
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(z,l+u/U)  endpoints  as  (0,1),  (-4.75,0.75)  and  (-8.5,1).  The 
area  subtended  by  this  approximation  and  the  vertical  axis 
1  +  u/U  -  1  gives 

0  *  1.06  cm  for  x  *  61.7  cm 

o 

In  Figure  (3b),  we  find  the  three  (z,l+u/U)  endpoints  as  (0, 
0.965),  (-4.5,  0.8)  and  (-9.75,  1).  The  subtended  area  gives 

0  ■  1.05  cm  for  x  ■  110.8  cm 

o 

In  Figure  (4),  the  velocity  profiles  predicted  by  equation 

(9)  using  the  estimated  value  0  -  1.05  are  plotted.  As  can  be 

o 

seen,  the  agreement  with  the  measured  profiles  is  Improved 
considerably  as  compared  to  Figure  3.  However,  the  predicted 
profiles  tend  to  be  shifted  downward  in  z  compared  to  the 
measured  ones  and  the  free  surface  velocities  are  somewhat  in 
error . 


The  Surface  Proximity  Factor 

If  we  assune  that,  because  of  the  proximity  of  the  free 
surface,  the  origin  of  uQ(z)  shifts  from  z  »  -d  to  z  -  -  ad, 
equation  (9)  becomes 
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He  term  a  the  surface  proximity  factor, 
plate  a  ■  1. 


For  a  deeply  towed 


To  evaluate  a  for  the  near  free  surface  experiment  of  Swean 
and  Keramidas  their  measured  free  surface  velocity  for  x  »  110.8 
cm  is  used.  In  particular,  referring  to  Figure  (4b),  we  require 
u(110.8,  0)/U  *  -0.05  which  provides  the  value  a  *  0.85. 

In  Figure  5,  the  velocity  profiles  predicted  by  equation 

(16)  with  0  -  1.05  and  a  »  0.85  are  shown.  The  agreement  with 

o 

the  measured  profiles  is  excellent. 


Conclusions 

Prandtl's  free  shear  layer  model  has  been  used  to  predict 
the  turbulent  wake  of  a  flat  plate  towed  below  a  free  surface. 
When  appropriate  account  is  taken  of  the  apparent  increase  in 
viscous  drag  on  the  plate  as  it  approaches  the  free  surface, 
excellent  agreement  between  predicted  and  measured  velocity 
profiles  is  obtained. 

These  results  validate  the  applicability  of  Prandtl's  model 
for  the  turbulent  kinematic  viscosity  for  this  complex  wake 
development  problem. 


Figure  3.  Free  shear  layer  prediction  from  eq .  (9)  of  the 

turbulent  wake  behind  a  flat  plate  towed  5  cm  below  a  free 

surface.  The  assumed  momentum  thickness  was  6  -  0.644  cm. 

o 

-  predicted,  ooo  measured  [I].  (a)  x  »  61.7  cm, 

(b)  x  ■  110.8  cm. 
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Figure  4.  Free  shear  layer  prediction  from  eq .  (9)  of  the 

turbulent  wake  behind  a  flat  plate  towed  5  cm  below  a  free 

surface*  The  assumed  momentum  thickness  was  6  *  1.05  cm. 

o 

-  predicted,  ooo  measured  (!].  (a)  x  *  61.7  cm. 


(b)  x 


110.8  cm 


1  +  u/U 


1  +  u/U 


Figure  S.  Free  shear  layer  prediction  from  eq 
turbulent  wake  behind  a  flat  plate  towed  5  cm 
surface.  The  assumed  momentum  thickness  was  6 
and  the  assumed  surface  proximity  factor  was  a 

-  predicted,  ooo  measured  [lj.  (a)  x  ■  61. 

(b)  x  ■  110.8  cm. 


.  ( 16 )  of  the 

below  a  free 

»  1.05  cm 
o 

-  0.85. 

7  cm , 


References 


1.  T.F.  Swean,  Jr.  and  G.A.  Keramidas  ,  1984:  The  Turbulent  Wake 
of  a  Flat  Plate  Near  a  Free  Surface.  Proceedings  of  the  2nd 
International  Conference  on  Computational  Methods  and 
Experimental  Mechanics,  Computational  Mechanics  Centre, 
Southampton,  England. 

2.  R.A.  Skop,  1984:  The  Hydrodynamic  Wake  of  a  Surface  Ship: 
Theoretical  Foundations.  Naval  Research  Laboratory  Report 
8833,  Washington,  DC.  AD  B084604L  DTT  June  28,  1984 

3.  R.A.  Skop,  1984:  The  Development  of  Plane  and  Initially 
Axisymmetrlc  Laminar  Wakes  Generated  by  Bodies  Moving  Beneath 
a  Free  Surface.  Proceedings  of  the  2nd  International 
Conference  on  Computational  Methods  and  Experimental 
Mechanics,  Computational  Mechanics  Centre,  Southampton, 
England . 

4.  H.  Schlichting,  1960:  Boundary  Layer  Theory,  Fourth  Edition. 
McGraw-Hill,  New  York,  NY. 

5.  B.A.  Boley  and  J.H.  Weiner,  1960:  Theory  of  Thermal  Stresses. 
John  Wiley  &  Sons,  New  York,  NY. 

6.  R.  Chevray  and  L.S.G.  Kovasznay,  1969:  Turbulence 
Measurements  in  the  Wake  of  a  Thin  Flat  Plate.  AIAA  Journal, 
Vol .  7,  No.  8,  pp.  1641-1643. 


